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We explore the possibility of the acoustic analogue of a super-radiance like phenomenon, i.e., 
the amplification of a sound wave by reflection from the ergo-region of a rotating acoustic black 
hole in the fluid draining bathtub model in the presence of the pressure be amplified or reduced 

in agreement with the value of the parameter ^7 = 1 + ^ . We remark that the interval of 

frequencies depend upon the neo-newtonian parameter 7 ( PLh = 777 LIh) and 

in this work. As a consequence, the tuning of the neo-newtonian parameter 

changes the rate of loss of the acoustic black hole mass. 


becomes narrow 



I. INTRODUCTION 

Acoustic black holes possess many of the fundamental properties of the black holes of general relativity and have been 
extensively studied in the literature [1-3]. The connection between black hole physics and the theory of supersonic 
acoustic flow was established in 1981 by Unruh [3] and has been developed to investigate the Hawking radiation 
and other phenomena for understanding quantum gravity. Hawking radiation is an important quantum effect of 
black hole physics. In 1974, Hawking combining Einstein’s General Relativity and Quantum Mechanics announced 
that classically a black hole does not radiate, but when we consider quantum effects emits thermal radiation at a 
temperature proportional to the horizon surface gravity. 

Since the Hawking radiation showed by Unruh [3] is a purely kinematic effect of quantum field theory, we can study 
the Hawking radiation process in completely different physical systems. For example, acoustic horizons are regions 
where a moving fluid exceeds the local sound speed through a spherical surface and possesses many of the properties 
associated with the event horizons of general relativity. In particular, the acoustic Hawking radiation when quantized 
appears as a flux of thermal phonons emitted from the horizon at a temperature proportional to its surface gravity. 
Many fluid systems have been investigated on a variety of analog models of acoustic black holes, including gravity 
wave [4] , water [5] , slow light [6] , optical fiber [7] and electromagnetic waveguide [8] . The models of superfluid helium 
II [9], atomic Bose-Einstein condensates [10, 11] and one-dimensional Fermi degenerate noninteracting gas [12] have 
been proposed to create an acoustic black hole geometry in the laboratory. 

The analogous systems employ a classical as well as Newtonian treatment generally, and also some quantum systems 
are considered. Also, a BH is a relativistic gravitational phenomenon which requires the inertial and gravitational 
effects of the pressure for a reasonable description of the system. However, relativistic pressure effects can be 
incorporated in a Newtonian framework in some cases approximatively. This is called neo-Newtonian theory and 
it is modification of the usual Newtonian theory by comprising the enough pressure into the dynamics. 

In the present work, our main goals is firstly to realize the process of drawing the acoustic BH for the neo-Newtonian 
hydrodynamics, on the other hand, to analyze the impact of neo-Newtonian parameter on the superresonance 
phenomenon especially on the frequencies of the waves. The paper is organized as follows. In Sec. II we will 
review the development process of the acoustic BH after recalling the basic concepts of the neo-Newtonian theory 
[13]. In Sec. Ill we will develop superresonance phenomenon and in Sec. IIIB the numerical results. In Sec. IV we 
make our final conclusions. 
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II. NEWTONIAN HYDRODYNAMICS IN AN EXPANDING BACKGROUND: COSMOLOGY 


In this section, we present the neo-Newtonian hydrodynamics applied to cosmology. However, we consider the 
standard case of Newtonian equations firstly. In the presence of inviscid perfect fluid, the basic equations of Newtonian 
hydrodynamics are 


0 , ( 1 ) 

-Vp, (2) 

where p, p and v are fluid density, corresponding pressure and the velocity field, respectively. The dot shows the 
differentiation with respect to cosmic time t. The above system of equations becomes suitable to study cosmology 
adopting the velocity field v = H(t)r (Hubble’s law) where H(t) = being a{t) the scale factor. It is worth noting 

the trivial solution for the continuity Eq. (1) p(a) = po/a 3 , where the today’s scale factor ao = 1 gives the today’s 
density of the fluid po- 

Gravitational interaction is coupled into Euler’s equation (2) as 

Vz? 

v + (v-V)v= - --VT, (3) 

P 

where the gravitational potential ’E obeys the Poisson equation 

V 2 ’E = 4t rGp. (4) 


p + V • {pv) = 
dv ■ 

p —j 7 = p[v+{v-V)v J = 


Eqs. (1) and (3) represents the fluid picture of cosmic medium which is gravitationally self-interacting through Poisson 
Eq. (4). In the framework of Newtonian cosmology, the Friedmann equations are given as follows 


(-2 E) 8t tG 


p and H + H 2 = — 


4t tG 


( 5 ) 


where E appears as integration constant associated to the energy of system. Moreover, the pressure does not 
correspond to homogeneous and isotropic background. With the inclusion of Newtonian cosmology, we can not 
model a radiation dominated as well as dark energy dominated epoch. This approach is restricted to a description of 
the Einstein-de Sitter universe. 


A. Including Pressure: The neo-Newtonian Theory 


The neo-Newtonian equations has been developed by McCrea [14] and Harrison [15] which ensures the effects of 
pressure as well the simplicity of Newtonian physics. Later, a crucial study related to the perturbative behavior of 
the neo-Newtonian equations has helped in developing the final expression for the fluid equations in this approach 
[16] (see also [17, 20]) which are given by 

d t pi + V • (piif) + pV • v = 0 . (6) 


VjD 

v + (v-V)v = -W -— , (7) 

P + P 

V 2 ^ = 47rG(p + 3p). (8) 

As is widely known, the above set of equations admits a homogeneous and isotropic solution, i.e. , p = pit) and 
p = p(t). In this case, the fluid velocity is (a dot means partial time derivative) 


a _ 
v = — r, 
a 


Combining Eqs. (9), (7) and (8), we can obtain the following equations 

a 2 (-2 E) _ 8irG 


H + H 2 = — ^^(p + 3 p). 


(9) 


( 10 ) 

( 11 ) 
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with the continuity equation ( 6 ) reducing to the following form 

^+3^(p + p) = 0. (12) 

These equations are exactly correspond to the relativistic Friedmann equations. The main idea behind the neo- 
Newtonian formalism relies on the following substitutions: Firstly, it is necessary to redefine the concept of inertial 
and passive gravitational mass density. With the redefinition 

Pi^p + p, (13) 


we rewrite the continuity and the Euler equation. 

The second step is the interpretation of the active gravitational mass density i.e., the density that source the 
gravitational field. Hence the following redefinition 


p g ^-p + 3p, (14) 

will become the source of the Poisson equation. The generalization of this result in the presence of pressure has been 
evaluated in [14]. Moreover, this approach has been modified [15] which leads to neo-Newtonian theory. 

When p = 0we find the Newtonian equations. The interesting feature of above equations is that the inertial mass 
present in the Newtonian equation can be replaced by p + From cosmological structure formation point of view, 
it is noted that equation ( 6 ) does not provide correct growth of matter for large scale perturbations in the scenario of 
homogeneous, isotropic and expanding background [18-20]. Actually, the correct growth for large scale perturbations 
is obtained if Eq.( 6 ) is modified as [18-20] 


d t p+ V • (pv) + pV • v = 0. (15) 

This equation shows correspondence with effective metric of usual Newtonian scenario when applied to fluid 
configuration considered. It represents that one can see, how the neo-Newtonian formulation is sensible to the 
the specific symmetries and hypothesis of the problem in another context of the cosmological one. This fact may 
indicate that the construction of a Newtonian counterpart of a relativistic problem may vary from problem to problem, 
deserving in the general case a deeper analysis. 


B. Acoustic Black Holes in neo-Newtonian Theory 


Let us now consider the barotropic fluid, i.e. p = p(p), inviscid and irrotational where the equation of state p = kp n , 
with k and n constants. We write the fluid velocity as v = — where ifj is the velocity potential. Thus, we linearise 
the equations ( 6 ), (7) and (15) by perturbing p, v and i/j as follows: 


P — Po + £/?i + 0(e 2 ) , 

p n = [po + epi + 0(e 2 )] ~ Pq + nepg 1 p\ + ... , 

v = Vo + ev\ + 0(er 2 ), 
i\) = i/j 0 + £ip + 0(e 2 ) , 

where p is the fluid density, p its pressure and v its the velocity field. Thus, the wave equation becomes 


- dt{c s 2 p 0 dt'tp + (- + ^)#o-VV> | + V • | - c s 2 p 0 v 0 (2 + 2 )^'Wo-'W’ + poV^j = 0, 
and we can rewrite as follows 

d^rd^) = 0 

where 




-1 

!+ 7 „ 


1 +7 7 .‘U 
2 u 


c s ~ (v x ) 2 — 7 v x v y 
-^p-v v —7 v x v v c 2 s — {v x ) 2 


(16) 

(17) 

(18) 
(19) 


( 20 ) 


( 21 ) 


(22) 
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Here we are assuming a (2+l)-dimensional spacetime. However, considering the Klein-Gordon equation [20] for a 
massless scalar field we obtain 


1 


= 0 , 


7=5 

and so the effective (acoustic) metric, in cylindrical coordinates, reads 
ds 2 = 


Po 


c 2 s + (v 2 +v 2 )(^rl 


— c 2 dt 2 + (dr — 7 v r dt) 2 (rd(j) + 7 dt) + dz z 


(23) 


(24) 


where 7 = 1 + —. 


C. Ergo-region and event horizon 

Now considering a static and position independent density, the velocity field is given by 


- A ~ B 1 

v = — r H- (p , 

r r 

which is a solution obtained from the continuity equation ( 6 ) and the velocity potential is 

ip(r, cf>) = — A In r — Bcp. 

Thus, considering (25), (26) and also the coordinate transformations as follows 

dt = dt' + a\dr, dip = dcj>' + a, 2 dr, 

being 

A r (1 + 7 ) A B 7 


a 1 = 


2 (A 2 7 — cs 2 r 2 ) ’ 2 r(A 2 7 — cs 2 r 2 ) ’ 

into the metric (24) we obtain the acoustic black hole in neo-Newtonian theory which is given by [13] 


ds 2 = poPi 


- 1 - 


dt 2 + (1 + fa) ( 1 — 


r 2 \ - 1 
r _h 
„2 


dr 2 - 


2 Bfo 


rdtdcp + ^1 + ^jr 2 d(p 


where 


A = (i + A)- 1/2 , A = (i (V) 


A = 


(1+7) 


04 = 


f A ( 7 - 1 ) 


being r e the radius of ergo-region and 77 the event horizon, i.e., 


r e = \/l( A2 + B 2 ), r h = ^\A\. 

Thus, the metric (29) can be now written in the form 

b/ 3 3 
r 

0 


g^iy — Pofdl 


-f 0 

0 (1 + 0 2 ) Q - 1 


B/3 3 


0 


(l+£). 


and the inverse of the g is 


<T = 


Pof3i(l + P 2 ) 


(1 + 7 t)Q _1 0 


0 

B/3 S 

rQ 


ZjlQ 

(i++F 


B03 

rQ 

0 

L 

Q 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 
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where 


/ 

~9 


= 1 -- 0 = 1 -^ 


(1 + fo) 
Q 


i + ^/ + ^ 


(35) 

(36) 


III. SUPERRESONANCE PHENOMENON 


A. Reflection coefficient 


Now, we introduce the metric (33) in the relation Klein Gordon given by 

^niV-Hg^d^ i) = o 


(37) 


yFa M 

Bearing in mind the stationarity and axisymetry spherical metric (33), we can separate the potential as follows: 

M*,r,0) =R(r) e i( ” f - ra#) (38) 

where m is a real constant (azimuthal number) and w the rotational frequency of the acoustic BH. Considering (33), 
(37), (38), The radial function R(r) satisfies the linear second order differential equation 


where 


R"(r) + Pi(r) R'(r) + P 2 (r) R(r) = 0 




(39) 


(40) 


and 


where 


and 


Ti = { C 2 r .2[ 8c 2 r 2 +jB 2 ( _ 1 + 7) 2 ]+ 3 A 4 7( _ 1+7 ) 2 +j 4 2 [ 3 B 2 ( _ 1+7)2 

X 1 + C 2 s ^(1 +7(6+7))]} 

T 2 = 2r[4c 2 r 2 + A 2 (—1 + q ) 2 + P 2 (—1 + 7 ) 2 ](c 2 r 2 — A 2 j) 2 


P 2 W = g 


r 3 = r 2 4c 2 + A 2 (-l + 7 ) 2 + B 2 (- 1 + 7) 2 


* 4 2 _2 1 2 2 

— 4c, r m + A r w 


(—1 + 7 ) + 4c 2 [(A 2 + B 2 )m 2 7 — i B m r u>(l + 7 ) + r 4 ttr], 


r 4 = 4 


4c 2 r 2 + A 2 (—1 + 7 ) 2 + P 2 (-l + 7 ) 2 ] (c 2 r 2 - A 2 7 ^. 


(41) 

(42) 

(43) 


(44) 


Thus, the problem has reduced to a one dimensional Schrodinger problem. Two further analytical simplifications 
can be made (39): first, We now introduce the tortoise coordinate r* by using the following equation 


A = (i_rni. 

dr* V r 2 ) dr' 


(45) 
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which gives the solution 


r 



(A 2 + B 2 ) 7 — 1 2 w \A\Sr 
4 c 2 1 2 ; ) \ 2 c s 


log 


r c s - A y/j 
r c s + A yfi 


)■ 


(46) 


Note that the tortoise coordinate spans the entire real line as opposed to r which spans only the half-line; the 
horizon r = ^ ^ maps to r* —> — oo, while r —> oo corresponds to r* —> -foo. Next, introducing a new radial 
function Z(r) G(r*) = R(r). The equation (39) becomes: 


d 2 G(r*) 

dr* 2 


+ Q(r) G(r*) = 0 


(47) 


with 


Q(r) = ^ ( 


1 ( 1 d 2 Z(r) , „ 1 dZ(r) , „ 

“l ' -;-h -T 2 


■Z’(r) dr 2 
(~l+(A 2 7 ]) 
(cs'rAy 


Z(r) dr 


) 


(144cs 4 r 2 ) 


16 ^ 1 + (a 2 +b 2 )(-i+ 7 ) 2 ^ 2 l (4cs 2 r 2 + A 2 (-l + y]) 2 + B 2 (-1 + y]) 2 ) 2 

+ 1__/2i8 CS v + d 2 r-i 

r 2 ( 4cs 2 r 2 + y4 2(_ 1+7 ])2 +j g2(_ 1+7 ])2) \ ^ ^ '• 

+ 7 ]) 2 + B 2 {- 1 + 7 ] ) 2 ) (cs 2 r 2 (8cs 2 r 2 + B 2 (- 1 + 7 ]) 2 ) + 3Al 4 (-l + 7 ]) 2 
x 7 ] + A 2 (3B 2 (-1 + 7 ]) 2 7 ] + cs 2 r 2 (l + 7 ] (6 + 7 ])))) j/ (r 2 (4cs 2 r 2 + 

x ,4 2 (-l + 7 ]) 2 + B 2 (- 1 + 7 ]) 2 ) 2 (cs 2 r 2 - ^ 2 7 ])| + (4r 2 (4cs 2 + A 2 
x (-1 + 7 ]) 2 + B 2 {- 1 + 7 ]) 2 )(-4cs 4 m 2 r 2 + A 2 r 2 {- 1 + 7 ]) 2 uf + 


4cs 2 ((A 2 + B 2 )m' 2 "f] — Bmr 2 ( 1 + 7])w] + r 4 w] 2 
x (-1 + 7]) 2 + B 2 (-1 + 7]) 2 )(csV - A 2 7]) 2 |. 


4cs 2 r 2 + A 2 


We remark that the main advantage of the new radial equation (47) over (39) is the absence in the former of a first 
derivative of the radial function. We analyze differential equation (47) in two distinct radial regions near the horizon, 
i.e r* —> —00 and at asymptote, i.e r* —> + 00 . In the asymptotic region, Eq.(47) can be written approximately as, 


d 2 G{r*) 
dr* 2 


+ ■ 


w* 


c 2 + l(A 2 + B 2 )(-l + 1 ) 2 


G(r*) = 0 


(48) 


of which one solution is 


G(r*) = exp 

+ 1Z exp 


w 


cf + j{A 2 + B 2 )(-l + ~/) 2 




— 1 r 


c 2 + ±(A 2 + B 2 )(-l + 7 ) 2 


(49) 


where 7Z is the reflection coefficient. In equation (49), the first term is the incident wave and the second term is the 
reflected wave. The Wronskian of solution (49) turns out to be 


W (+ 00 ) = —2 i i 


c l + t(A 2 + B 2 ){-1 + 7) 2 


( 1 -|^| 2 ). 


(50) 


Considering the second field (r* —> — 00 ), the equation (47) becomes: 

d 2 G{r*) 


dr* 2 


+ {4(-2 B c 2 m + ^ 2 (1 + 7) w) 2 }/{A 2 (4c 2 r 2 + A 2 (-1 + 7 ) 2 































+ B 2 {-1 + 7 ) 2 )(H 2 (1 + 7 ) 2 + B 2 (-1 + 7 ) 2 )}G(r*) = 0. 
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The physical solution of equation (51) is the following. 

G(r*) = T exp — ir*|4(— 2.B c 2 m + A 2 (l + 7 ) w ) 2 } j jdl 2 (4c 2 r 2 

+ A 2 (- 1 + 7 ) 2 + B 2 (- 1 + 7) 2 )(A 2 (1 + 7 ) 2 + S 2 (-l + 7 ) 2 )}. 


where T" is the transmission coefficient. The Wronskian of solution (52) is following 


W{— 00 ) = 

0 

+ S 2 (-l + 7 ) 


— 2 z 4(—21? c 2 m + A 2 (1 + 7 ) w) 

,2 ^2 1 / i 2 / 1 1 _.\2 1 n 2 / 1 1 ..\2 




x (4cf r 2 + H 2 (-l + 7) 2 + B 2 {-1 + 7 ) 2 )(A 2 (1 + 7 ) 2 


(51) 


(52) 


(53) 


Since both equations are actually limiting approximations of the differential Eq. (47), which, as we have mentioned, 
has a constant Wronskian, it follows that ([21], [22], [23], [24]) 

W(— 00 ) = W(+oo), (54) 

so that, from (50) and (53), we obtain the relation 

r 1 1 25 m c s 1 

\TZ\ 2 = 1- , - =T 2 . (55) 

^/(A 2 (l+ 7 ) 2 + B 2 (-l + 7 ) 2 ) 

superresonance phenomenon occurs when the norm of the reflected wave is greater than the norm of the incident 
wave, that is, when the reflection coefficient is greater than unity [25, 26]. Hence, we can observe in eq. (55) that, for 
frequencies in the range 


0 < w < m fin 


(56) 


with 


Ah 


2 

1 + 7 


Ah, 


Ah = 



(57) 


the reflection coefficient has a magnitude larger than unity whose imply the amplification relation of the ingoing sound 
wave near horizon, with this condition we can extract the energy of the system [24]. It is noteworthy that 7 > 1 
makes the interval to have superradiance smaller. Here m is the azimuthal mode number and Ah is the angular 
velocity of the usual Kerr-like acoustic BH. the angular velocity of the usual Kerr-like acoustic BH depends of the 
pressure (parameters n and k). Thus, we show that the presence of the the neo-newtonian parameter 7 modifies the 
quantity of removed energy of the acoustic BH and that is either possible to accentuate or attenuate the amplification 
of the removed energy of the acoustic. The effect of superresonance can be eliminated when 7 = (— 1 ). 


B. Numerical results 


In order to confirm the effect of pressure on the phenomenon of superresonance, it is necessary to extract the 
potential energy. The objective in this section, would be to see if the pressure changes the potential energy In this 
way we obtain the following differential equation for the radial function R{r) (39) as 


0 + §) 


2 2Bf3^rrnjj m 2 f 

r * j w ^2 ^2~ 


(1 + j3 2 )R{r) 

(—5)2 


We can rewrite the equation (58) as 




u} 2 — 


2Bf3^mijj m 2 f 


R(r) 


1 d 
r^J—g dr 

J-(r) d 

” l - 

r dr 


rV^g( l + /?2) 1 Q- r 


dr 


dr 


R(r) = 0 , 


R(r) = 0. 


(58) 


(59) 

























where T(r) = y/—g{ 1 + P 2 ) 1 Q(r). At this point we introduce the coordinated g using the following equation [27] 


d_ 

dg 



(60) 


and now introducing the new radial function G(g) = r 1 / 2 i?(r), we can obtain the following modified radial equation 
we get to a new radial equation obtained from (59) that reads 


d 2 G(g) 

dg 



Bfem 

r 2 




G{g) = 0 , 


(61) 


and the potential V (r) is given by [29, 30] 


V 0 ) 


F( r ) 4 m 2 f(r) _ 

4r 2 F(r) 1 ’ 


i -E^( 1 + ?±y l + 2r d -Sd 

r A J~ (r) \ r z J dr 


(62) 


We note that the equation (62) does not satisfy the asymptotic behavior V(r) —> 0 as r —> 00 . We can see from 
Fig.(l) when the 7 parameter becomes larger, the potential energy increases, which allows us to say that the presence 
of pressure ascents the phenomenon of superradiance. 



2 4 6 8 10 


y=l ,9(Neo-Newtonian) 


y= 1,6(Neo-Newtonian) 


y= 1.3(Neo-Newtonian) 


y=l (Newtonian) 



FIG. 1: The graphs illustrating the evolution of the potential V(r) for each value of 7 . The function is plotted for 
A = B = c s =m = 1. 


One important aspect of the results reported here is that, perhaps, the specific form of the neo-Newtonian equations, 
retaining the essential features of the relativistic problem, may depend on the symmetries of the problem, as a 
comparison with the cosmological case suggests. We have constructed acoustic black hole configurations in the 
context of the neo-Newtonian hydrodynamics. The presence of the parameter 7 in the effective metric of acoustic 
BHs in the framework of the neo-Newtonian theory is the main feature of our results. This parameter is always larger 
than 1 if k and n have the same sign, that is, if the fluid have the expected positive square speed of sound. Hence, 
essentially, the velocity on the horizon is diminished by the presence of the pressure. It is worth noting that when 7 
is equal to 1, the results of Newtonian theory are recovered [1, 3], which leads us to rewrite 7 as follows 


lneonewtonian 'Ynewtonian 


kn Pq 1 


(63) 


The supplementary term is the contribution of the pressure in neo-Newtonian theory. 




















9 


IV. CONCLUSION 

In this paper we shown that the presence of the pressure modify the quantity of removed energy of the acoustic 
BH and that, it is possible to accentuate or to attenuate the amplification of the removed energy of the acoustic BH 
and still exists the possibility to cancel the superradiance effect i.e the reflec- tion coefficient is equal to unity, when 
7 = ( 2t "^ 1h — i) where SI# is the angular velocity of the acoustic black hole. Furthermore, the interval of frequencies 
becomes narrower in the present scenario. As a consequence, the tuning of the neo-newtonian parameter 7 changes 
the rate of loss of mass of the acoustic BH. 

Acknowledgement: Authors thank Prof. J. C. Fabris and Prof. H. E. S. Velten for useful comments and 
discussions. Ines G. Salako thanks IMSP for hospitality during the elaboration of this work. 
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